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Liquid crystals are phases of matter intermediate between crystals and liquids. Whereas classical 
liquid crystals have been known for a long time and are used in electro-optical displays, much less is 
known about their quantum counterparts. There is growing evidence that quantum liquid crystals 
play a central role in many electron systems including high temperature superconductors, but a 
quantitative understanding is lacking due to disorder and other complications. Here, we analyse 
the quantnm phase diagram of a two-dimensional dipolar gas, which exhibits stripe, nematic and 
supersolid phases. We calculate the stiffness constants determining the stability of the nematic 
and stripe phases, and the melting of the stripes set by the proliferation of topological defects is 
analysed microscopically. Our results for the critical temperatures of these phases demonstrate 
that a controlled study of the interplay between quantum liquid and superfluid phases is within 
experimental reach for the first time, using dipolar gases. 


The investigation of cold atomic gases has enabled one 
to study many-body physics with unrivalled experimen¬ 
tal control and in regimes never realised before. Recent 
progress in trapping and cooling of gases consisting of 
dipolar atoms/molecules opens up a promising new re¬ 
search direction. The dipole-dipole interaction is long 
range and anisotropic, which is predicted to give rise to a 
number of exotic forms of matter [1, 2]. Degenerate Fermi 
gases consisting of atoms with a large magnetic dipole 
moment have already been created [3, 4], and progress to¬ 
wards producing degenerate gases of fermionic molecules 
with an electric dipole moment is being reported [5, 6]. 

Here, we analyse the quantum phases of a two- 
dimensional (2D) dipolar Fermi gas at non-zero temper¬ 
atures T. This includes a stripe phase, whose low en¬ 
ergy degrees of freedom are described by an anisotropic 
XY model. We determine the stiffness constants of this 
effective model microscopically. The finite temperature 
melting of the stripe phase is driven by the proliferation 
of topological defects called dislocations, and the corre¬ 
sponding Berezinskii-Kosterlitz-Thouless (BKT) critical 
temperature is determined by the well-known renormal¬ 
isation group equations. For large tilting angles of the 
dipoles, the system can have additional superfluid pair¬ 
ing which coexists with the stripe order. We calculate the 
critical temperature of the superfluid transition. When 
the dipoles are perpendicular to the 2D plane, the critical 
temperature of the stripe phase is shown to vanish, and 
the system exhibits a nematic phase characterised by long 
range orientational order but no translational order. Our 
results demonstrate that with dipolar gases, it is within 
experimental reach to study quantum liquid and super¬ 
fluid phases characterised by varying degrees of spon¬ 
taneous translational, rotational, and gauge symmetry 
breaking. The interplay between such phases is believed 
to play an important role in many electronic materials 
discovered in recent decades [7-11]. Moreover, our results 
show that one can confirm the microscopic mechanism 


behind the BKT transition, namely the proliferation of 
topological defects, simply by observing the proliferation 
of disclocation defects in the stripe pattern. Such an ex¬ 
perimental verification of BKT physics has been achieved 
only recently using atomic gases [12], whereas other ex¬ 
periments reported only indirect evidence of BKT physics 
in the bulk properties [13-20]. 


RESULTS 

We consider fermionic dipoles of mass m and aver¬ 
age areal density uq, which are restricted to move in 
the xy plane by a tight harmonic trapping potential 
nvJlz^ 12 along the z-direction. In the limit ep, 

where ep = kp/2m = 2TTno/m is the Fermi energy of 
a 2D non-interacting gas with areal density no, the sys¬ 
tem is effectively 2D with the dipoles frozen in the har¬ 
monic oscillator ground state in the z direction. An ex¬ 
ternal field aligns the dipoles so that their dipole mo¬ 
ment d is perpendicular to the y-axis and forms an an¬ 
gle 0 with the z-axis. The dipole-dipole interaction is 
Vd(r, z) = D^(l —3cos^ (^rd)/, where is the 
angle between the relative displacement vector of the two 
dipoles (r, z) with r = (x, y) and the dipole moment d, 
and = d?fAneo for electric dipoles and = (Pyo/Air 
for magnetic ones. 

The strength of the interaction is determined by the di¬ 
mensionless parameter g = AmD^kp and the degree 
of anisotropy is controlled by the tilting angle 0. The 
system is rotationally symmetric for 0 = 0 and becomes 
more anisotropic with increasing 0. Above a critical in¬ 
teraction strength gd^), it is predicted to form density 
stripes at T = 0, where the density exhibits periodic 
modulations of the form 

n(r) = no + ni cos(qc ■ r — u). (1) 
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Here, qc || is the wave vector of the stripes, and 
ni and u are their amplitude and phase respectively. 
The density modulation is formed along the j/-direction 
so as to minimise the interaction energy. The system 
thus exhibits liquid-like correlations along the x-direction 
and crystalline correlations along the y-direction. This 
phase has been predicted by Hartree-Fock theory [21- 
25], density-functional theory [26], and by a variant of 
the so-called STLS method [27]. Remarkably, Hartree- 
Fock and density-functional theory predict essentially the 
same critical coupling strength gc(0) ~ 0.6 for stripe for¬ 
mation at 0 = 0, whereas the STLS method obtains a 
somewhat higher value. For 0 > 0.237r, the system is 
predicted to become a p-wave superfluid [28], which for 
strong enough coupling can coexist with the stripe order 
forming a supersolid [37]. The dipoles are also predicted 
to form a Wigner crystal for y ~ 27 for 0 = 0 [29, 30]. 
This very strong coupling regime is outside the scope of 
the present paper. 

Stripe phase at finite T and effective XY model 

Since the stripe phase breaks translational invariance 
along the y-direction, it is a quantum analog of a classical 
smectic liquid crystal [31, 32]. Indeed, the system has a 
manifold of equivalent ground states distinguished only 
by a constant factor u, which specifies the position of the 
stripes along the y-direction. Consequently, there are low 
energy collective excitations associated with a spatially 
dependent phase u{x,y). Moreover, since a change from 
M to u -I- 27 r returns the system to the same ground state, 
it follows that the low energy degrees of freedom of the 
stripe phase are described by a 2D anisotropic XY model. 
Specifically, the simplest form of the elastic free energy 
congruent with the symmetry of the system is given by 

Fei = lJ + Bi^idyuf] = I y d\{Wuf 

( 2 ) 

for 0 yf 0. Here, B± and i?|| are the perpendicular and 
parallel elastic coefficients describing respectively the en¬ 
ergy cost of small rotations and compressions/expansions 
the stripes. In the second equality, we have used the 
rescaling x —>■ y/B\\/B±x to obtain an isotropic XY 
model with the effective elastic constant B = 

Berezinskii-Kosterlitz-Thouless melting 

As the stripe phase is described by the XY model, 
it exhibits algebraic long-range order at sufficiently low 
temperatures and it melts via the Berezinskii-Kosterlitz- 
Thouless mechanism due to the proliferation of topolog¬ 
ical defects [33-36]. In the case of the stripe phase, the 
topological defects are dislocations. The phase field for 


a single dislocation of charge Q = ±I,±2... satisfies 
f Vu(r) ■ dl = 27r(5, where the path of the integration en¬ 
closes the core of the dislocation. The presence of such a 
dislocation corresponds to inserting Q extra stripes to the 
left (right) of the dislocation for Q > 0 (Q < 0). The en¬ 
ergy of a single dislocation consists of a core part Ec, and 
a part that scales logarithmically with the size of the sys¬ 
tem. Pairs of bound dislocations with opposite charges 
(5 = ±I (IQl > I are energetically suppressed), however, 
have a finite energy even for an infinite system size and 
can be thermally excited in the stripe phase. This is 
due to the fact that the phase fields of the oppositely- 
charged dislocations cancel at large distances, which re¬ 
sults in merely a local disturbance of the density. In Fig. 
I, we illustrate dislocation pairs with opposite charges 
Q = ±I centered at (xi,yi) and (x 2 ,j/ 2 ) respectively. 
The stripe amplitude is suppressed in the core regions of 
the defects due to the large energy cost associated with 
Vt6 oc I/r, where r is the distance to the core. From the 
rescaling x —>■ ,yB^\/Bj_x it follows that the energy of 
a vertically displaced dislocation pair distance 6 apart is 
the same as that of a pair displaced horizontally by the 
distance ^yB±/B^^S. Since Ry > B± as we will demon¬ 
strate below, this shows that the dislocation pairs along 
the x-direction are more tightly bound than those along 
the j/-direction. 


The spontaneous thermal excitation of bound disloca¬ 
tion pairs decreases the elastic coefficients at a macro¬ 
scopic scale. The softening of the effective stiffness con¬ 
stant B can be calculated from the well-known renor¬ 
malisation group equations as described in the methods 
section. At a critical temperature T®*, the renormalised 
elastic coefficient drops to zero by a sudden jump of 
magnitude 2T®*/7r. This disappearance of elastic rigidity 
signals the melting of the density stripes. 

Calculation of bare stiffness constants 

We now turn to a microscopic calculation of the “bare” 
stiffness constants B\\ and B± unrenormalised by disloca¬ 
tion pairs. The relevant thermodynamic quantity is the 
free energy of the system F{q), which depends on the 
stripe wave vector q. Any non-uniform phase fluctuation 
increases the free energy by an amount given by (2) for 
long wave lengths. To extract the elastic coefficients Bj_ 
and i3||, we consider two specific distortions: an infinitesi¬ 
mal rotation and an infinitesimal compression/expansion 
of the stripes away from the equilibrium configuration, as 
illustrated in Fig. 2. 

These distortions are described by the phase field u(r) = 
—(5q • r, where (5q = (5(7|| • By for the compression and 
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Figure 1. Dislocation pairs in the stripe phase. The dis¬ 
locations are centered at {xi,yi) and (x 2 , 2 / 2 ) so that the phase 
field is M(r) = arctan[(y — t/i)/(a; — xi)] — arctan[(j/— i/2)/(x — 
X 2 )]. Top: A Q = 1 dislocation centered at (0, —2.1-e jqc) and 
a Q = —1 dislocation centered at (0, 2.77r/5c)- Bottom: A 
Q = 1 dislocation centered at {—Sn/qc, 0) and a Q = —1 dis¬ 
location centered at (37r/gc,0). The dashed lines indicate the 
position of the density maxima. 


Figure 2. Elementary distortions of the stripes. Rota¬ 
tion (left) and compression (right) of the stripes away from 
their equilibrium positions indicated by the solid lines. 


i5q = dq±-ex for the rotation. They are thus equivalent to 
a variation of the stripe vector q = qc-l-<5q. Inserting the 
phase fluctuations into Eq. (2), we obtain the increment 
of the free energy 

Fei = F(q, + dq) - E(q,) = (3) 


where A is the area of the system, and we have used the 
equilibrium condition VF(q)|q^ = 0. We thus find 




1 d^F{q) 
A dql 


Bn 


1 d^F(q) 
A ~d^ 


( 4 ) 


The interaction energy per particle due to stripe forma¬ 
tion scales as ^^^(ni/nq)^. Assuming that the in¬ 

teraction energy is dominant, we find that the elastic 
coefficient B scales as ~ (ni/no)^(/eF for a fixed 0. The 
magnitude of B can be further reduced by a geometrical 
factor depending on 0 , since the system becomes rota- 
tionally symmetric for 0 = 0 , as we shall discuss below. 

In order to microscopically calculate the bare stiffness 
constants, we employ Hartree-Fock mean-field theory for 
the free energy, writing F(q) ~ TMF(q) = 
where Hmf is the mean-field thermodynamic potential 
given by 


f^MF = -Fint - T^ln [1 + 


( 5 ) 


/i is the chemical potential and JV is the total number 
of particles. The quasiparticle energies are where 
j = 1, 2 , • • • is the band index and k is restricted to the 
first Brillouin zone of the ID periodic potential set up by 
the stripes. We subtract the interaction energy Eint to 
avoid double counting. The details of this calculation are 
given in the methods section. 

In Fig. 3, we plot the bare elastic coefficients obtained 
from this approach as a function of temperature for 5=1 
and 0 = 0 . 287 r, for which the system has a large stripe 
amplitude ni/no ~ 0.6 at low temperatures. In order 
to minimize finite size effects, we determine the elastic 
coefficients by fitting a parabolic curve to the free energy 
in the vicinity of q = qc in accordance with (3), instead 
of performing a numerical differentiation following (4). 
This is illustrated in the insets of Fig. 3. This procedure 
allows us to obtain numerically accurate values for the 
elastic coefficients. From Fig. 3, we see that both elastic 
coefficients decrease with increased temperature. This is 
expected since thermal excitations of quasi-particles re¬ 
duce the stripe amplitude and thus their rigidity. We 
also find that B± By, which suggests that compress¬ 
ing/expanding the stripes costs more energy than a ro¬ 
tation. This difference in magnitude becomes even more 
profound for small 0 when B± is strongly suppressed by 
the weak anisotropy of the system. Finally we note that 
for 5 = 1 and 0 = 0 . 287 r, the system is in fact predicted 
to have additional superfluid pairing at T = 0 [37]. How¬ 
ever, as demonstrated in Ref. [37], the superfluid order 
has negligible effects on the stripe formation, and it can 
thus be safely neglected when analysing the elastic prop¬ 
erties of the stripes. 
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Figure 3. Bare stiffness constants as a function of tem¬ 
perature. Upper curve gives and lower curve -Bj_ for 
g = 1 and 0 = 0.287r. The blue circles in the insets are elas¬ 
tic free energy plotted as a function of q at T = O-OIcf, and 
the red solid curves are parabolic fits to several data points 
in the vicinity of q = qc. 

In Fig. 4, we plot the bare elastic constants as a func¬ 
tion of the tilting angle 0 for 5 = 1 and T = O.Olei?- 
The elastic constant i?|| depends non-monotonically on 
0 , first decreasing and then increasing exhibiting a min¬ 
imum at 0 ~ 0.247r. This is consistent with the mean- 
field phase diagram, which shows that the stripe forma¬ 
tion is somewhat suppressed for intermediate values of 
0 [25, 37]. To illustrate this, we plot as an inset the 
stripe amplitude ni as a function of 0 ; we see that it ex¬ 
hibits the same non-monotonic behaviour as By. In com¬ 
parison to this behaviour. Fig. 4 shows that B± increases 
monotonically in 0. In particular, we have B± —>■ 0 for 
0 —)■ 0 as shown in detail in the inset. This reflects that 
the system is rotationally symmetric for 0 = 0 such that 
a rotation of the stripes costs no energy. 


Renormalised stiffness constants and stripe melting 

The bare elastic constants obtained from the mean- 
field theory can now be used as initial values in the RG 
equations to determined the renormalised elastic con¬ 
stants. We also need the dislocation core energy, which 
must scale as ~ R. Therefore, we write E^. = kB, where 
K is a constant of order unity. In Fig. 5, we plot the 
renormalised elastic coefficient Bn as a function of tem¬ 
perature, obtained by solving (9) with the initial mean- 
field values ot B — -^JBYB^i and Ec = kB for various 
coupling strengths g and tilting angles 0. To examine 
the dependence on the core energy, we have chosen dif¬ 
ferent values of k. We see that the thermal excitation 
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Figure 4. Bare stiffness constants as a function of tilt¬ 
ing angle. Upper curve gives i3|| and lower curve B± for 
g — 1 and T — O.OIef. For B±, the red dashed line is an ex¬ 
trapolation for 0 < O.IStt, where the coefficient is too small 
to be accurately determined with our numerical method. The 
upper inset is a plot of the relative stripe amplitude as a func¬ 
tion of 0. The lower inset is an expanded view of the B± for 
small values of 0. 

of dislocation pairs soften the elastic coefficients as ex¬ 
pected. This softening is negligible for low T where the 
core energy prohibits the excitation of dislocations. The 
softening increases with decreasing core energy and in¬ 
creasing T. At the critical temperature T®* determined 
by the solution to 2 r/ 7 r = Bn{T), the elastic coefficient 
drops to zero discontinuously and the stripes melt. 



Figure 5. Renormalised stiffness constants as a func¬ 
tion of temperature. The three groups of curves for BjiiT) 
for g = 1 correspond to, in the order from bottom to up, 
0 = 0.27r, 0 = 0.267r and 0 = 0.287r respectively. In each of 
the group, the four curves, in the order from bottom to up, 
correspond to k = 2, 3, 4 and 00 (mean-field result) respec¬ 
tively. The slope of the dashed diagonal line is 2/7r. 
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The resulting melting temperature T®* is plotted in 
Fig. 6 as a function of 0 for g = 1 and k = 3. It increases 
rapidly with 0, indicating that the degree of anisotropy 
of the system increases such that the stripes become more 
rigid. An extrapolation of our calculation for g = 1 and 
K = 3 shows that T®* ^ O.OGci? for 0 ~ 0.37r. The critical 
temperature also increases with the coupling strength, 
scaling as Tf ^ B {ni/no)'^gep. We note that in 
addition to the explicit linear dependence on g, the T®* 
can further increase with the coupling strength through 
the dependence on ni. Our results show that in order to 
observe the stripe phase and the associated BKT physics 
with dipoles, it is preferable to choose a large tilting angle 
in addition to having a large dipole moment. However, 
the tilting angle cannot exceed 0 ~ O.Stt above which 
the system exhibits a density collapse for large coupling 
strengths [28, 37]. 



Figure 6. Phase diagram for g — 1. The system is in the 
stripe phase below the Kosterlitz-Thouless melting temper¬ 
ature which is calculated taking k = 3 for the core energy. 
For 0 = 0, the striped phase melts at T = 0+ into a nematic 
phase with long range orientational order but no translational 
order. The nematic phase melts into an isotropic liquid at 
a temperature T" ~ indicated by the blue cross. For 
0 > arcsin(2/3) the system is in a supersolid phase at T = 0 
with both stripe and superfluid order. The transition temper¬ 
atures of this phase calculated for various tilting angles are 
indicated by the diamonds. The dashed curve is a fit to the 
data by (7). 


Melting of supersolid phase 

The system exhibits p-wave pairing for 0 > 

arcsin(2/3) [28], which can coexist with stripe order for 
g > <7c(0) at r = 0 [37]. We now determine the critical 


temperature Tf for the superfluid transition. The 2D su¬ 
perfluid transition is in principle also determined by the 
BKT mechanism, where the topological defects are now 
vortices. For weak pairing, however, the mean-field BCS 
theory in fact gives a good estimate of the transition tem¬ 
perature. We thus determine the critical temperature by 
solving the linearised gap equation 

/2T) r 
2Ck' 

( 6 ) 

Here Ak is the gap parameter and V(k, —k') is the effec¬ 
tive interaction between the quasiparticles in the stripe 
phase with energy dispersion ^k measured from the Fermi 
surface. The details of this calculation are given in the 
methods section. The critical temperature obtained from 
this calculation is shown in Fig. 6 for 5 = 1 and for sev¬ 
eral tilting angles. This mean-field result gives an upper 
bound to the critical temperature, but since Tf <C ef we 
expect that a more detailed BKT calculation yield only 
slightly smaller values. This should be contrasted with 
the melting of the stripes, where an estimate of the criti¬ 
cal temperature from a vanishing stripe order would give 
a much higher value compared to the BKT calculation. 
This can be seen from Fig. 3, which shows that the mean- 
field elastic coefficients remain large up to T = 0.126^. 
Thus, it is crucial to use the BKT theory to analyse the 
stripe melting. 

Using a simple p-wave ansatz for the gap parameter 
Ak — A cos (f), where (j) is the polar angle of the wave 
vector k, one can obtain an approximate solution for the 
critical temperature as 

Tf (7) 

where C is a constant related to an effective momentum 
cutoff in the integral in (6). We find that the data ob¬ 
tained from solving (6) numerically are in fact very well 
described by (7) with C ~ 0.4. 

Quantum nematic phase for 0 = 0 

Figure 6 shows that the critical temperature for the 
stripe phase vanishes as 0 —)■ 0. This is a direct conse¬ 
quence of the rotational symmetry rendering B±^ = 0 for 
0 = 0. In this case, the system is no longer described 
by the XY model. Instead, an appropriate expression for 
the elastic energy of stripe fluctuations is [38] 

Fe\ = J (frlX'^idluf + {dyuf], ( 8 ) 

where A is a length comparable to the stripe spacing. 
Dislocations again play an important role in determining 
the finite temperature properties of the system described 


Ak = — 


dk' 

(27r)' 


■V(k,-k')Ak- 


tanh(^k' 


2Ck' 
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by (8). In contrast to the 6 > 0 case, however, single 
dislocations now have a finite energy and can be ther¬ 
mally excited. When the presence of the free dislocations 
is taken into account, a system described by (8) is pre¬ 
dicted to be in a nematic phase for 0 < T < T", and in 
an isotropic liquid phase for T > T" [38]. In the nematic 
phase, the translational order exists only within a length 
scale which is determined by the density of the free 
dislocations. The stripe orientations, averaged over the 
length scale are however algebraically correlated. As 
a crude physical picture, one can think of the nematic 
phase as blobs of stripe order of area ~ , which are all 

oriented more or less in the same direction, but which are 
not positionally correlated with each other. The nematic 
phase is in this sense analogous to the 2D hexatic phase 
of a crystal, which exhibits bond orientational order but 
no long-range translational order [32, 39, 40]. A quan¬ 
tum hexatic phase was recently predicted to exist in 2D 
dipolar gases for very strong coupling g ^ 27 [41, 42]. 
The results presented here point out the intriguing pos¬ 
sibility to realise a quantum version of the nematic phase 
with dipoles for smaller coupling strengths. We expect 
the critical temperature T" for the melting of the quan¬ 
tum nematic phase to scale as By. However, a quantita¬ 
tive calculation of the critical temperature for the dipolar 
system requires knowledge of the parameter A, whose de¬ 
termination is beyond our current theoretical framework. 
In Fig. 6, we have indicated the critical temperature T” 
using a somewhat smaller value than the bare Hy due to 
renormalisation effects. 


DISCUSSION 

An important question concerns whether the critical 
temperature for the predicted quantum liquid crystal 
phases is within experimental reach. As an example, let 
us consider a recent experiment reporting the trapping 
of chemically stable ^^Na^°K molecules in their ground 
state close to quantum degeneracy. The group obtained 
an induced dipole moment of d = 0.8Debye and a maxi¬ 
mum 3D density of usd = 2.5 x 10^^cm“^ [6]. Estimat¬ 
ing a corresponding 2D areal density as uq = , these 

values correspond to g ~ 0.57. This coupling strength 
can be increased by reaching a larger fraction of the 
permanent electric dipole moment of ^^Na‘*°K, which is 
2.72Debye [43], or by increasing the density of the gas. 
Since the critical temperature for the nematic and the 
stripe phases both scale as ~ {ni/noYgeF, this indicates 
that the quantum liquid crystal physics discussed in this 
paper is within experimental reach, once dipolar gases 
can be cooled down significantly below their Fermi tem¬ 
perature. 

The formation of stripe and superfluid order can be 
observed as correlation peaks in time-of-flight (TOF) ex¬ 
periments [25, 37]. One can also detect the stripes di¬ 


rectly as density modulations, either after TOF or in-situ, 
provided that the experimental resolution is sufficiently 
high. Observing the proliferation of dislocations would 
directly confirm the microscopic mechanism behind the 
BKT transition. 

Finally, we would like to mention a recent fixed note 
Monte-Carlo calculation which suggests that the striped 
phase is not the ground state for 0 = 0 for any coupling 
strength [44]. We speculate that this result, which con¬ 
tradicts that of Refs. [21-27], is due to the approximate 
nature of the calculation combined with the fragility of 
the striped phase, which melts at any non-zero tempera¬ 
ture for 0 = 0, as shown by our results. 

In summary, we analysed the phase diagram of a 2D 
dipolar gases, which exhibits stripe, nematic and super¬ 
solid phases corresponding to the breaking of transla¬ 
tional, rotational and gauge symmetry. For a non-zero 
tilting angle 0, the low energy degrees of freedom of 
the striped phase are described by an anisotropic 2D 
XY model. We calculated the stiffness constants cor¬ 
responding to a rotation and a compression/expansion 
of the stripes microscopically. This should be contrasted 
with electron systems, where such stiffness constants are 
often simply unknown parameters of the theory. The 
stripes were shown to melt via the Berezinskii-Kosterlitz- 
Thouless mechanism due to the proliferation of disloca¬ 
tions, and we obtained the melting temperature using the 
relevant renormalisation group equations. We also calcu¬ 
lated the critical temperature of the supersolid phase. 
For 0 = 0, the striped phase is stable only at T = 0, 
which melts into a nematic phase for arbitrarily small 
temperatures. Our analysis of the melting temperatures 
demonstrated that they should be within experimental 
reach. An observation of these phases would constitute 
a major breakthrough in our understanding of the inter¬ 
play between liquid crystal and superfiuid order in low¬ 
dimensional many-body systems. 


METHODS 

Renormalisation group equations 

We calculate the softening of the effective stiffness con¬ 
stant B = due to the excitations of dislocation 

pairs using the well-known renormalisation group equa¬ 
tions 


dK ^(l) _ D s„,2/ 


dl 


= 27rYil) 


dyjl) 

dl 


= [2-7rK{l)]y{l). 


(9) 


Here K{1) = B{l)/T and y{l) = exp[—ifc(l)/T] are the 
scale-dependent stiffness constant and dislocation fugac- 
ity respectively. They both decrease with increasing I 
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as the renormalisation due to dislocation pairs at larger 
length scales are included via the solution of (9). The 
initial values of K{0) and y{0) are the bare (local) values 
unrenormalised by dislocation pairs, which we calculate 
microscopically as described in the text. At a critical 
temperature T®*, the long range renormalised elastic co¬ 
efficient Bji = limi_>oo 7?(?) drops to zero by a sudden 
jump of 2r®‘/7r. This disappearance of elastic rigidity 
signals the melting of the stripes. 


Mean-fleld theory of stripe formation 


The mean-field Hamiltonian that takes into account 
the possibility of stripe formation with a wave vector q 
is given by [25] 

iiMF = ^ ekcj^Ck + y^[/ikCk+a4 + h.c.], (10) 

k k 

where c|^ creates a dipole with momentum k, Ck is the 
single particle Hartree-Fock energy 

ek = ^ - k')](c[,Ck-), (11) 

k' 

and /^k is a real ofF-diagonal element defined by 

^ - k')](c|,,Ck-+q). (12) 

k' 

The quasi-2D interaction in Fourier space is obtained 
by averaging the interaction over the harmonic oscilla¬ 
tor ground state in the 2 direction. This gives (up to an 
irrelevant constant term) [45] 

V (k) ~ ^ (cos^ 0 — sin^ 0 cos^ (p) , (13) 

where (p is the polar angle of k. We diagonalise the mean- 
field Hamiltonian by generalising the method described 
in Refs. [25, 37] to an arbitrary stripe vector q. This 
yields the Hamiltonian 

BmF ^jk'ljk^jk* 

Here y^k = Fj k+cCk+c annihilates a quasiparticle 
with energy where j = 1,2, ••• is the band index, 
G = ?q, I = 0, ±1, • • • is the reciprocal lattice vector and 
k is restricted to the first Brillouin zone of the ID periodic 
potential set up by the stripes. We can then calculate the 
mean-field free energy as FMF(q) = + P-N, where 

riiviF is the mean-field thermodynamic potential given by 
(5) and N = fjk with = [exp/3(£^k - A^) + 1]"^- 


The interaction energy is most easily calculated using 

^'MF = ^jkfjk — ^kin + 2i?i„t (15) 

i.k 


where Akin = X]k(^k^k)fc^/2TO the kinetic energy. 


BCS theory of the superfluid transition 


To explore superfluid pairing within the stripe phase, 
we use BCS theory with the quasiparticle Hamiltonian 
’Hbcs = TAmf + 'Hp. Here, 


Hp = 


E Evk'^]'.-k')%-kE-k + h.c. 

ii'kk' 


describes pairing between the time-reversed quasiparti¬ 
cles, interacting via 


Vy,(k',-k) 


E '^G-G'.G'-cE.k'+CJ'E.-k'+G' 

GG'GG' 


xE-k+cEk+G^ (k _ k' + G - G') . (16) 


To derive a gap equation that is amenable to a par¬ 
tial wave expansion, we switch to the “extended zone 
scheme”, whereby a single particle state ipjkip) in the 
j’th band in the first BZ is mapped onto a state 4’kip) 
in the j’th BZ in the standard way [37], where the vector 
k is now unrestricted. The effective pairing interaction 
Vj'j{k',—k) shall be denoted by V(k, —k') and quasi¬ 
particle dispersion by Sk- Pairing between time- 
reversed quasiparticles gives rise to the gap parameter 
Ak = X^k'h’(k, —k')(7_k'7k'), which satisfies the finite 
temperature gap equation 


/2T) V 
2Ck' 
(17) 

Here ^k = £k — P and Ak = + l^kp, where the 

chemical potential p is approximated by the value in the 
stripe phase. The Cauchy principal value term 'P/2^k' in 
(17) renders the gap equation well defined with no need 
for a high momentum cut-off. At temperatures in the 
vicinity of the superfluid transition, the linearisation of 
the above gap equation yields (6) we use in the main text. 
Equation (6) can be solved by the method of partial wave 
expansion described in Ref. [37]. Finally we determine 
the transition temperature by gradually increasing T in 
the gap equation until it ceases to admit finite solutions. 


Ak = — 


dk! 

{2f)‘ 


■V(k,-k')Ak 


tanh(Ak 


2Ak 
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